Astronomy & Astrophysics manuscript no. vega-vl4 


©ESO 2013 


January 20, 2013 





Letter to the Editor 



High-sensitivity search for clumps in the Vega Kuiper-belt 

New PdBI 1.3 mm observations ** 

V. Pietu 1 , E. Di Folco 2 - 3 , S. Guilloteau 4 5 , F. Gueth 1 , and P. Cox 1 



o 

(N 



(N 



in 

6 



> 

oo 
in 

in 
o 



X 



IRAM, 300 rue de la piscine, F-38406 Saint Martin d'Heres, France e-mail: pietu@iram. fr 

CEA Saclay/Service d'Astrophysique, Laboratoire AIM, CEA/DSM/IRFU-CNRS-Universite Paris Diderot, F-91191 Gif-sur- 
Yvette Cedex, France 

LESIA-UMR 8109, CNRS and Observatoire de Paris-Meudon, 5 place J. Janssen, 92195 Meudon, France 

Universite de Bordeaux, Observatoire Aquitain des Sciences de l'Univers, 2 me de l'Observatoire BP 89, F-33271 Floirac, France 
CNRS/INSU - UMR5804, Laboratoire d'Astrophysique de Bordeaux; 2 rue de l'Observatoire BP 89, F-33271 Floirac, France 



Received ; accepted 



ABSTRACT 



Context. Previous studies have found that Vega is surrounded by an extended debris disc that is very smooth in the far infrared, but 

displays possible dumpiness at 850 jim and dust emission peaks at 1.3 mm. 

Aims. We reobserved Vega at 1.3 mm with PdBI to constrain its circumstellar dust distribution. 

Methods. Our observations of a three-field mosaic have a factor of two higher sensitivity than previous observations. 

Results. We detect Vega photosphere with the expected flux, but none of the previously reported emission peaks that should have 

been detected at the > 6<r level, with a sensitivity <lrhK. 

Conclusions. This implies that the dust distribution around Vega is principally smooth and circularly symmetric. This also means 
that no planet is needed to account for dust trapped in mean-motion resonnance. 



Key words. Radio continuum: stars - Stars: Vega - Stars: circumstellar matter - Method: observational 
eters 



Instrumentation: interferom- 



1. Introduction 

More than 25 years a fter the first debris discs were discov ered 
by the IRAS mission dAumann etalj|1984t lAumann|[T985h . the 
presence of Edgeworth-Kuiper belts around other stars has been 
widely attested. Statistical studies conducted with the Spitzer 
satellite on large samples of main sequence stars have shown 
that 30 + 5 % of A type and 16 + 3 % of FGK ty pe stars exhibit 
detectable IR excess emission longward of 24yum dSu et alj2 006; 
iBrvden et alj |2006: Trill ing et al.l i2008). This excess luminosity 
is attributed to the presence of cold dust grains arranged in a disc 
or ring(s). Among those, only a handful of nearby systems have 
been resolved because of obvious limitations in sensitivity and 
angular resolution. Being one of the first systems detected, Vega 
(A0V, 7.76 pc) has a particular status, and is usually considered 
as a prototype of debris discs stars. 

Given its proximity to the Earth, its face-on disc has been 
successfully imaged from IR to mm- wavelengths. Its appear- 
ance curiously seems to change from a smooth azimuthal pro- 
file at short wavelengths to a more structured, clumpy ri ng in the 
sub-m m domain. MIPS images at 24, 70 and 160/im dSu et alj 
120051) show circular profiles, the emission being detected up 
to 815 AU for the latter. The radial profile analysis suggests a 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 

** Processed data from Fig[T]are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



ring-like distribution of material (after subtraction of the pho- 
tosphere contribution), with an inner boundary located around 
85 ± 15 AU. Herschel-PACS observations at 70 and 160/mi di- 
rectly confirmed t his ring morphology, with a peak intensity oc- 
curing at 85 AU (ISibthorpe et alJGOlOl) . as well as the smooth 
brightness profile. 



At sub-mm wav elengths, the first image was reported by 
iHolland et all dl998l) with JCMT/SCUBA at 850/mi. The re- 
solved emission extents to scales twice as large as our Solar 
System, with a possible peak emission located at ~ 70 AU, 
northeast of the star, accounting for less than 40 % of the to- 
tal flux. The possible asymmetry of their map, although at the 
2<r level, has motivated several subsequent studies to explore 
the nature of the blob a nd its putative link with unseen plan- 
ets. iKoerner etall feoOl) observed Vega with OVRO at 1.3 mm 
and reported a ring arc located at 95 AU (12"from the star). 
IWilner et al.l d2002l) used the Plateau de Bure Interferometer 
(IRAM) to probe the continuum emission at 1.3 and 3.3 mm. 
While the star was only detected at the former wavelength, 
no peak emission appeared in the origi nal resolutio n map . 
Degrading the interferometer resolution, Wil ner et al.l (P002) 
claimed that extended exists emission concentrated in two blobs 
located at 8 and 9.5"(60 and 75 AU), southwest and northeast 
of the star r e specti vely. Using the SHARC II camera at CSO, 
iMarsh et"a"D d2006) detected ring-like structures of the 350/vm 
and 450/mi emissions, with some inhomogeneity at the 2-4 cr 
level. The brightness peaks in the photosphere-subtracted im- 
ages partially overlap with those reported at 850//m and 1.3 mm, 



2 



V. Pietu et al.: High-sensitivity search for clumps in the Vega Kuiper-belt 



Field 


RAOff. 


Dec Off. 


Noise 


Fit. RA 


Fit. Dec 


Flux 


Corr. Flux 




(") 


(") 


(mJy) 


(") 


(") 


(mJy) 


(mJy) 


North-East 


3 


5 


0.26 


0.13± 0.16 


0.30±0.14 


1.70+0.23 


2.05+0.28 


South-West 


-3 


-5 


0.29 


-0.12± 0.19 


0.03±0.16 


1.48±0.23 


1.81+0.28 


Central 








0.37 


-0.70± 0.15 


-0.01±0.16 


1.96+0.32 


1.98±0.32 



Table 1. Description of the fields observed and point-source fit. From left to right, the columns are: Col. 1: Field name. Col.2,3: 
Field pointing center (offsets from a =18:36:56.330, 5 =38:47:01.30, J2000.0). Col. 4: Noise in the field. Col. 5,6,7,8: Fit of point- 
source. Col. 5,6: Fitted positions are offsets from a =18:36:56.359, 5 =38:47:01.20 (after correction for proper motion - see text). 
Col. 7,8: Fitted flux, and flux corrected for primary beam attenuation. 



strongly suggesting that grains and/or planetesimal concentra- 
tions exist in the belt around Vega. 

The possible existence of concentrations in the dust belt has 
long been interpreted (and modelled) as the evidence of grains 
trapped in mean-motion resonances (MMR) with a putative, but 
unseen, massive planet located close to the inner radius of the 
plane tesimal ring (iHolland et al.lll998LlWilner et al.l2002tlWvatil 
120031 Two scenarios have been invoked to populate the MMRs 
with an unseen planet: the particles are trapped either while 
migrating inwar d under the influence of P oyinting-Robertson 
(P-R) drag (e.g. iKuchner & Holmanl 12003). or during th e out- 
ward m igr ation of the p l anet i tself. Models developed by Wvatt 
(2003) and lReche et alJ (12008b invoke the outward migration of 
a Neptune- to Saturn-mass planet in a circular (and elliptical, 
respectively) orbit to reproduce the observed clumps. To recon- 
cile the different views of the disc morphology at the different 
wavelengths, it has been suggested that the clumps seen in the 
sub-mm trace large grains trapped in the same resonance as their 
parent planetesimals, while intermediate and small grains (down 
to the blown-out size) are decoupled from the resonant struc- 
ture, thus forming a smooth and axisymmetric disc component 
observed at shorter wavelengths. This broad and smooth halo 
could therefore be sculpted by small gr ains evacuate d from the 
MMRs by the stellar radiation pressure (Wvatt 2006) and/or by 
the high-velocity c ollisions between their parent planetesimals 
dKrivov et alj |2007). Substantial observational efforts have been 
invested in the quest for the predicted planet(s), mostly based 
on near-IR adaptive optics (and/or coronagraphic) tech niques 
on large telescopes (Ma rois et al.ll2006t iHinz et al.l l2006). This 
search remains unsuccessful today, essentially because of the 
lack of sensitivity at the predicted planet position (7 - 10"). 

The Vega system thus appears to be a unique laboratory for 
constraining and understanding the evolution of dust grains in 
exo-Kuiper belts, in particular the interaction between planetes- 
imals, grains, and planets. 



2. Observations and analysis 

2.1. Observations 

During the winter 2006/07, a new generation of receivers was 
installed on the Plateau de Bure interferometer, enabling single 
frequency dual-polarization observations with improved noise 
figures, single-side band tuning and larger bandwidth. This in- 
strumental upgrade resulted in a significantly improved sensitiv- 
ity of the array (at 1 .3 mm, the total gain was an increase by a fac- 
tor of four in continuum sensitivity). As part of the science veri- 
fication program, Vega was observed in January 2007 at 230.538 
GHz. These observations were carried in the D (compact) con- 
figuration of the PdBI (with baselines ranging from 15 to 85 m, 
the most suited arrangement for detecting large scale emission). 
The tuning was in the lower side-band (LSB), with T rec ~ 50 K 



and a bandwith of 850 MHz in each of the two orthogonal lin- 
ear polarizations (H and V). The atmospheric phase noise ranged 
from 20 to 50°. Despite these first observa t ions r eaching a sen- 
sitivity comparable to that of IWilner et al.l (120021) . we detected 
only the central point source, but no extended emission at all. 

We thus reobserved a two-field mosaic towards Vega (offset 
by (3 , 5") and (-3 , -5"), respectively, the primary beams being 
indicated in Fig Q] together with the central pointing obtained in 
January) in October and November 2007, again in the D con- 
figuration and with the same spectral setup. Atmospheric phase 
noise was again between 20° and 50°. The absolute flux cali- 
bration was performed using observation of MWC349. Its flux 
density model has recently been revised at IRAM, with a 16% 
increase in flux (Krips et al. in prep), so we used 1.98 Jy as a 
flux model at 230.5 GHz. With a grand-total observing time of 
16 hours on source, we finally reach a sensitivity limit below 
0.2 mJy/beam (0.9 mK with natural weighting) in the central re- 
gion of the final mosaic map, and < 0.4mJy/beam within 11" 
from Vega. 

2.2. Analysis 

Vega has a high proper motion (jj = 201.03 mas/yr, fi s = 
287. 47 mas/yr, with an erro r ellipse of 0.63 x 0.54 mas/yr, at PA 
144°. lPerrvman et al.ll997h . We correcte d the individual (u, v) ta- 
bles for proper motion using the value of Perrv man et al.l(11997l) 
to a common J2000.0 Equinox prior to data analysis. 

We then performed point-source fits in the (u,v) plane and 
detected a point source in each of the fields. The fitted position 
agrees within the error bars (see Table [I], as does the flux den- 
sity, when we take into account the attenuation by the primary 
beam for the offset fields. The mean value is 1.95 + 0.17 mJy. 
IWilner et al.l (12002 ) found 1.7+0.30mJy ft but used 1.70 Jy as 
a flux density for MWC349. Correcting for the calibrator flux 
difference, their flux becomes 1 .98+0.35 mJy, in excellent agree- 
ment with our measurement. This is consistent with the 2.11 mJy 
value expected fro m the pole-on photosp here for which we refer 
to the modelling by Mii ller et al.l (1201 Ol) . who took into account 
the discrepancy between the polar and equatorial temperature re- 
sulting from the fast spin velocity of this star. 

Data were processed with GILD AS softwarj^] using the clas - 
sical mosaicing algorithm (see e.g. iGueth & Guil loteau 2000). 
We present the resulting continuum maps in Fig. [TJ in the top 
row we present maps obtained with natural weighting, while in 
the bottom row the data have been tapered in the (u,v) plane 
(with elliptical Gaussian 1 je widths of 30 x 40m at PA 0°) to in- 
crease the sensitivity to putative extended emission by degrading 

1 They cite 1.7 + 0.13 mJy from a Gaussian fit, but the error in a 
point source flux cannot be below the rms noise of the image, which 
they indicate is 0.3 mJy/beam 

2 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about GILD AS. 
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Fig. 1. Top left: natural weighting 1.3 mm continuum image (corrected for the primary beam attenuation). The angular resoluion 
is 2.5x2.1" at PA 82° and the contour spacing 0.5mJy/beam (2.1 mK), with negative contours being dashed. The dashed circle 
represents the primary-beam half -power field of view. Top center: corresponding noise map. Contour spacing is 0.2mJy/beam . 
The circles show the observed individual fields (half-power field of view, centered on the crosses). Top right: signal-to-noise ratio 
map. Contour spacing is 2 cr. Bottom left: 1.3 mm continuum image obtained with a (u, v) taper. Angular resolution is 4.2x3.7" at 
PA 60° and contour spacing is 0.6 mJy/beam (0.8 mK). Bottom center: corresponding noise map. Contour spacing is 0.2 mJy/beam 
(with the first contour being at 0.4mJy/beam). Bottom right: signal-to-noise ratio map. Contour spacing is 2cr. In all maps, the 
filled square symbols indicate the location o f the extended emissi ons found bv lWilner et al.l (120021) and the empty square symbols 
indicate the location of the peaks reported bv lKoerner et al.l (1200 lh . North is up, and east to the left. 



the angular resolution. Since the mosaicing algorithm produces 
images corrected for the primary beam attenuation, the noise in- 
creases towards the map's edges (as is clearly visible in the sig- 
nal maps in the left column of Fig. [TJ. We present in the middle 
column of Fig. [T]the corresponding noise maps that show our 
sensitivity limit as a function of the position in the field of view. 
The right column of Fig. Q] shows the signal-to-noise ratio maps 
obtained in each case by dividing of the signal map by the noise 
map. From these, we conclude that only the central point source 
is detected with a SNR> 3. 

2.3. Comparison with previous work 

Pr evious interferometri c ob servations at 1.3 mm were reported 
bv lKoerner eTail (1200 ll) and lWilner etall d2 002). With an angu- 
lar resolution of 3.3 x 2.9", Koe rner et al.l d2001l) found an arc 
ring at 95 AU with four peaks located between 11 and 14.5" 
of flux density ranging from 2.4 to 4.0 mJy (3 to Act level) af- 
ter correcting for pr imary beam attenua tion. In addition to the 
central point source. IWilner et all (120021) reported two emission 
peaks located at 8.0 and 9.5", after tapering their data to a 5.3 
x 4.6" resolution. From a Gaussian fit in the (u,v) plane, they 
obtained fluxes (corrected for the primary beam attenuation) of 



7.1 + 1.4 and 4.3 + 1.0 mJy, respectively (4 to 5<x level). The 
position of these dust peaks is reported in Fig. Q~| with their lcr 
error ellipse derived from the signal-to-noise ratio and synthe- 
sized beam. These two results clearly do not agree, as none of 
the error ellipse overlap. We also note that the comparison of the 
different results need to take into account that their single field 
maps are not corrected for primary beam attenuation, in contrast 
to the maps we present here. 

At the position of the peaks reported by iKoerner et al.l 
d200lh . our observations have a lcr sensitivity of 0.35 mJy/beam 
(0.54 for the most distant peak), a factor of between two and 
three better than the OVRO observations. If real, these peaks 
should have been detected at the > 6 - 8<r level. Similarly, 
our observations reach a lcr sensitivity limit of 0.26 mJy/beam 
(1.1 mK) or 0.45m Jy/beam (0.6 mK) at the position of the 
IWilner et al.l (|2002) peaks for the untapered and tapered maps, 
respectively. Thus, depending on their size, the two structures 
should have been detected at the 10 - 15<x and 6 - 10cr level in 
our maps. 

Among the possible ways of explaining our non-detection, 
we exclude calibration errors since we detect the central source, 
at the expected position and with the expected flux, which is 
an excellent internal quality check. Given the above-mentioned 
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flux calibration scaling with respect to the IWilner et al.l ([2002) 
data, any supposed bias in the flux calibrationshould enhance 
the probability of detecting the blobs in our maps, since their 
intensity should also be higher by 16%. If these structures had 
rotated according to Keplerian motion, the displacement within 
the ~ 6 year interval would have been smaller than 1". 



3. Discussion 

Our PdBI observations are consistent with no detectable milli- 
metric emission from the outer dusty disc (or ring) surrounding 
Vega. We briefly discuss the impact of this new constraint on the 
dust properties in this system. 



3.1. dumpiness 

We have shown in section I2.3I that t he clum ps p reviously re- 
ported at 1.3 mm by Koe mer et al.l (1200 ll) and IWilner et al.l 
(2002) are not detected in our observations, which are deeper in 
sensitivity by a factor of two, and are probably low-significance 
artifacts. A trivial implication of this finding is that no planet 
is required to trap dust in mean-motion resonance at these 
positions. More generally, our results also shed new light on 
the dumpiness of the dust surrounding Vega. Observations at 
shorter wavelengths display very smooth images dSu et alj 2005; 
Sibth orpe et al.ll20Toh and the asymmetry of the 850yum image 
is only at the 2< r level, hence also c ompatible with a smooth 
dust distribution (Holla nd et al.l[l99 8). The only detection of a 
possible asy mmetry comes fro m the SHARC II images at 350 
and 450 ^m dMarsh et al.ll2006h but is at the 2-4 ex and < 10% of 
the ring brightness. In conclusion, all observations display or are 
compatible with a smooth dust distribution. 



3.2. Upper limit on the disc emission at 1.3 mm 

Earlier sub-mm observations have convincingly established the 
existence of an extended emission around Vega, with a ring- 
like morpho l ogy a t 350/zm, 450 jum with CSO SHARC II 



(Mar sh et all I2006I). a nd 70/mi-160jUm with Herschel/PACS 



Sibthorpe et al. 2010I) . According to these studies, the peak in- 



tensity is located around 85 - 100 AU (1 1 - 13"). To set an upper 
limit on the 1.3 mm emission of this ring, we fitted a thin, uni- 
form, 85 AU radius ring into the (u,v) data, carefully taking into 
account the different fields of the mosaic (for numerical reasons, 
this thin ring is represented by an adequately sampled series of 
point sources). We find an integrated flux of 0.1 + 2.0 mJy, i.e. 
a 3 cr upper limit of 6 mJy. The lack of clear detection at high 
angular resolution relative to the bolometer results, which are 
more sensitive to extended flux, suggests that a thin ring may 
not be the most appropriate distribution. If we instead assume a 
uniform annulus of inner radius 85 AU and outer radius 105 AU, 
we find an integrated flux of -5.3 +2.4 mJy, i.e., a 3cr upper limit 
of 7.2 mJy, which is similar to the narrow ring value. 

Using the Sibtho rpe et all d2010h model, with an inner 
Gaussian and outer exponential decline in surface brightness 
(with a transition radius of 14"), our results are compatible with 
a total flux < 24 mJy (3 cr) if the distribution is truncated at 200 
AU, and < 30mJy if it extends up to 300 AU. This modelling 
also agrees with the 3cr upper limit of 6 mJy for emission within 
105 AU. With this model, we are clearly unable to place strong 
constraints on the dust properties, since our data lack sensitivity 
to faint extended emission. 



Assuming a dust temperature of 68 K (ISheret et al.l2004l) an d 
a dust absorption coefficient of 0.5 cm -g 1 (iNatta et alT f2004). 
our 3 cr flux density limit of 6 mJy translates into an upper limit 
on the dust grain mass in the ring of Md us t = 0.8M moon (or 
10~ 2 Mffl, not t aking into account the parent bodies). With the 
Sibthorpe et al.l ( 120101) model, the upper limit on the total dust 
grain mass is 3.9 M moon . This upp er limit is consisten t with the 
maximum dust mass obtained by iMiiller et al.l d2010l) in their 
modelling of the Vega disc as a steady-state collisional cascade. 

The best fit of t he dust spectra l energ y di stribution us- 
ing m easurements of Sibthor pe et al.l ([2010) and iHolland et al.l 
(Il998h . corrected for the stellar emission, gives a spectral index 
of 2.9 and results in an extrapolated flux at 1.3 mm of 12mJy. 
Our upper limit hence agrees with the measurements made at 
shorter wavelengths. 



4. Conclusions 

Taking advantage of new receivers at IRAM Plateau de Bure 
Interferometer, we have probed the Vega debris disc to investi- 
gate in depth the possible presence of clumps in the well-known 
dust belt. Apart from the stellar photosphere, we have detected 
neither a point-source nor extended emission in our three-field 
mosaic, although our sensitivity limit at 1.3 mm has improved 
by a factor of two compared to earlier studies. In particu- 
lar, we have found no evide nce for the nort hea st and south- 
west b lobs claimed by Koe rner et alj <2001l) or IWilner et alj 
(2002). These structures, which had been interpreted as the sig- 
nature of the gravitational influence of putative unseen planets, 
should have been detected in our interferometric map at the > 6cr 
level. 

The upper limits that we have derived on both the emission 
of a narrow ring extended emission are compatible with obser- 
vations at shorter wavelenghts. All observations indicate that the 
disc/ring system surrounding Vega is smooth and circularly sy- 
metric, hence do not require the gravitational influence of a giant 
planet to sculpt asymmetries in the disc, at the current level of 
detection. 

Acknowledgements. We acknowledge help from the Plateau de Bure staff when 
carrying out the observations. We thank our anonymous referee for helpful com- 
ments. 
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